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Abstract. We consider the system of (> 2) 
elastically colliding hard balls of masses mi, . . . ,mN 
and radius r on the flat unit torus T'^, u > 2. We 
prove the so called Boltzmann-Sinai Ergodic Hypoth- 
esis, i. e. the full hyperbolicity and ergodicity of 
such systems for every selection (mi, . . . , m^v; of 
the external parameters, provided that almost every 
singular orbit is geometrically hyperbolic (sufficient), 
i. e. the so called Chernov-Sinai Ansatz is true. The 
present proof does not use the formerly developed, 
rather involved algebraic techniques, instead it em- 
ploys exclusively dynamical methods and tools from 
geometric analysis. 
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§1. Introduction 

In this paper we prove the Boltzmann-Sinai Ergodic Hypothesis under the con- 
dition of the Chernov-Sinai Ansatz (see §2). In a loose form, as attributed to L. 
Boltzmann back in the 1880's, this hypothesis asserts that gases of hard balls are 
ergodic. In a precise form, which is due to Ya. G. Sinai in 1963 [Sin(1963)], it states 



1 Research supported by the National Science Foundation, grants DMS-0457168, DMS-0800538. 

1 



2 



that the gas of > 2 identical hard balls (of "not too big" radius) on a torus T^, 
v >2 (a zv-dimensional box with periodic boundary conditions), is ergodic, provided 
that certain necessary reductions have been made. The latter means that one fixes 
the total energy, sets the total momentum to zero, and restricts the center of mass 
to a certain discrete lattice within the torus. The assumption of a not too big radius 
is necessary to have the interior of the configuration space connected. 

Sinai himself pioneered rigorous mathematical studies of hard ball gases by prov- 
ing the hyperbolicity and ergodicity for the case N = 2 and u = 2 in his seminal 
paper [Sin(1970)], where he laid down the foundations of the modern theory of 
chaotic billiards. Then Chernov and Sinai extended this result to (A^ = 2, z/ > 2), 
as well as proved a general theorem on "local" ergodicity applicable to systems of 
N > 2 balls [S-Ch(1987)]; the latter became instrumental in the subsequent stud- 
ies. The case N > 2 is substantially more difficult than that of N = 2 because, 
while the system of two balls reduces to a billiard with strictly convex (spherical) 
boundary, which guarantees strong hyperbolicity, the gases of > 2 balls reduce to 
billiards with convex, but not strictly convex, boundary (the latter is a finite union 
of cylinders) - and those are characterized by very weak hyperbolicity. 

Further development has been due mostly to A. Kramli, D. Szasz, and the present 
author. We proved hyperbolicity and ergodicity for N = 3 balls in any dimension 
[K-S-Sz(1991)] by exploiting the "local" ergodic theorem of Chernov and Sinai [S- 
Ch(1987)], and carefully analyzing all possible degeneracies in the dynamics to 
obtain "global" ergodicity. We extended our results to = 4 balls in dimension 
u > 3 next year [K-S-Sz(1992)], and then I proved the ergodicity whenever N <v 
[Sim(1992)-I-II] (this covers systems with an arbitrary number of balls, but only in 
spaces of high enough dimension, which is a restrictive condition). At this point, 
the existing methods could no longer handle any new cases, because the analysis 
of the degeneracies became overly complicated. It was clear that further progress 
should involve novel ideas. 

A breakthrough was made by Szasz and myself, when we used the methods of 
algebraic geometry [S-Sz(1999)]. We assumed that the balls had arbitrary masses 
mi, . . . ,mN (but the same radius r). Now by taking the limit — >^ 0, we were 
able to reduce the dynamics of N balls to the motion of A^ — 1 balls, thus utilizing a 
natural induction on A^. Then algebro-geometric methods allowed us to effectively 
analyze all possible degeneracies, but only for typical (generic) (A^ + l)-tuples of 
"external" parameters (mi, . . . , rriN, r); the latter needed to avoid some exceptional 
submanifolds of codimcnsion one, which remained unknown. This approach led to 
a proof of full hyperbolicity (but not yet ergodicity!) for all A^ > 2 and > 2, and 
for generic (mi, . . . , m^v, r), see [S-Sz(1999)]. Later the present author simplified 
the arguments and made them more "dynamical" , which allowed me to obtain full 
hyperbolicity for hard balls with any set of external parameters (mi, . . . ,mN,r) 
[Sim(2002)]. (The reason why the masses rrii are considered geometric parameters 
is that they determine the relevant Riemannian metric 
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of the system, see §2 below.) Thus, the hyperbohcity has been fuUy estabhshed for 
all systems of hard balls on tori. 

To upgrade the full hyperbohcity to ergodicity, one needs to refine the analy- 
sis of the aforementioned degeneracies. For hyperbohcity, it was enough that the 
degeneracies made a subset of codimension > 1 in the phase space. For ergodic- 
ity, one has to show that its codimension is > 2, or to find some other ways to 
prove that the (possibly) arising codimension-one manifolds of non-sufficiency are 
incapable of separating distinct ergodic components. The latter approach will be 
pursued in this paper. In the paper [Sim(2003)] I took the first step in the direction 
of proving that the codimension of exceptional manifolds is at least two: I proved 
that the systems of iV > 2 disks on a 2D torus (i.e., i/ — 2) are ergodic for typical 
(generic) (A^ + l)-tuples of external parameters (mi, . . . , mjv, r). The proof again 
involves some algebro-geometric techniques, thus the result is restricted to generic 
parameters (mi, . . . jUIn; r). But there was a good reason to believe that systems 
m u > 3 dimensions would be somewhat easier to handle, at least that was indeed 
the case in early studies. 

Finally, in my paper [Sim(2004)] I was able to further improve the algebro- 
geometric methods of [S-Sz(1999)], and proved that for any N > 2, u > 2 and for 
almost every selection (mi,...,mAr; r) of the external geometric parameters the 
corresponding system of N hard balls on T'^ is (fully hyperbolic and) ergodic. 

In this paper I will prove the following result. 

Theorem. For any integer values N > 2, u > 2, and for every {N + l)-tuple 
(mi, . . . ^tunit) of the external geometric parameters the standard hard ball sys- 
tem {yiffi^ri I'S'^r}' A*m,r^ is (fuUy hyperbolic and) ergodic, provided that the 
Chernov-Sinai Ansatz (see §2) holds true for all such systems. 

Remark 1.1. The novelty of the theorem (as compared to the result in [Sim(2004)]) 
is that it applies to each {N -\- l)-tuple of external parameters (provided that the 
interior of the phase space is connected), without an exceptional zero-measure set. 
Somehow the most annoying shortcoming of several earlier results has been exactly 
the fact that those results are only valid for hard sphere systems apart from an 
undescribed, countable collection of smooth, proper submanifolds of the parameter 
space ]R-^+-^ 3 (mi,m2, . . . jUIn; r). Furthermore, those proofs do not provide any 
effective means to check if a given (mi, . . . , m^v; r)-system is ergodic or not, most 
notably for the case of equal masses in Sinai's classical formulation of the problem. 

Remark 1.2. The present result speaks about exactly the same models as the 
result of [Sim(2002)], but the statement of this new theorem is obviously stronger 
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than that of the theorem in [Sim(2002)]: It has been known for a fong time that, for 
the family of semi-dispersive biUiards, ergodicity cannot be obtained without also 
proving full hyperbolicity. 

Remark 1.3. As it follows from the results of [C-H(1996)] and [0-W(1998)], all 
standard hard ball systems (M, {S'^jt^K, //) (the models covered by the theorem) 
are not only ergodic, but they enjoy the Bernoulli mixing property, as long as they 
are known to be mixing. However, even the K-mixing property of semi-dispersive 
billiard systems follows from their ergodicity, as the classical results of Sinai in 
[Sin(1968)], [Sin(1970)], and [Sin(1979)] show. 

The Organization of the Paper. In the subsequent section we overview the 
necessary technical prerequisites of the proof, along with many of the needed ref- 
erences to the literature. The fundamental objects of this paper are the so called 
"exceptional manifolds" or "separating manifolds" J: they are codimension-one 
submanifolds of the phase space that are separating distinct, open ergodic compo- 
nents of the billiard flow. 

In §3 we prove our Main Lemma (3.5), which states, roughly speaking, the fol- 
lowing: Every separating manifold J C M contains at least one sufficient (or geo- 
metrically hyperbolic, see §2) phase point. The existence of such a sufficient phase 
point a: G J, however, contradicts the Theorem on Local Ergodicity of Chernov and 
Sinai (Theorem 5 in [S-Ch(1987)]), for an open neighborhood U of x would then 
belong to a single ergodic component, thus violating the assumption that J is a 
separating manifold. In §4 this result will be exploited to carry out an inductive 
proof of the (hyperbolic) ergodicity of every hard ball system, provided that the 
Chernov-Sinai Ansatz (see §2) holds true for all hard ball systems. 

In what follows, we make an attempt to briefly outline the key ideas of the proof 
of Main Lemma 3.5. Of course, this outline will lack the majority of the nitty-gritty 
details, technicalities, that constitute an integral part of the proof. 

We consider the one-sided, tubular neighborhoods Us of J with radius (thickness) 
5 > 0. Throughout the whole proof of the main lemma the asymptotics of the 
measures n{Xs) of certain (dynamically deflned) sets Xs C Us are studied, as 
5 — > 0. We fix a large constant C3 3> 1, and for typical points y E Us\Us/2 (having 
non-singular forward orbits and returning to the layer Us\Us/2 infinitely many times 
in the future) we define the arc-length parametrized curves Py,t{s) (0 < s < h{y,t)) 
in the following way: py^t emanates from y and it is the curve inside the manifold 
Sq(|/) with the steepest descent towards the separating manifold J. Here Sq^/) is 
the inverse image (E^(y)) of the flat, local orthogonal manifold (flat wave front, 
see §2) passing through y^ — S^{y). The terminal point H(j/) = ^ {h{y,t)) of the 
smooth curve t is either 

(a) on the separating manifold J, or 
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(b) on a singularity of order ki = ki{y). 
The case (b) is further split in two sub-cases, as follows: 
(b/1) ki{y)<cs; 
(b/2) C3 < ki{y) < oo. 

The set of (typical) points y e Us \Us/2 with property (a) (this is the set 
Us{oo) in §3) is handled by lemmas 3.28 and 3.29, where it is shown that, actually, 
Us{oo) = 0. Roughly speaking, the reason of this is the following: For a point 
y G Us{oo) the powers of the flow exhibit arbitrarily large contractions on the 
curves py^t (Appendix II), thus the infinitely many returns of S^{y) to the layer 
Us \ Us/2 would "pull up" the other endpoints (n.{y)) to the region Us \ J, 
consisting entirely of sufiicient points, and showing that the point n(y) e J itself is 
sufficient. 

The set Us\Us{cs) of all phase points y G Us\Us/2 with the property ki{y) < C3 
are dealt with by Lemma 3.27, where it is shown that 

fi{Us\Us{c3))=o{6), 

as 5 — > 0. The reason, in rough terms, is that such phase points must lie at the 
distance < 5 from the compact singularity set 

U s-'{sii-), 

0<t<2c3 

and this compact singularity set is transversal to J, thus ensuring the measure 
estimate p {Us \ Usics)) — o{5). 

Finally, the set Fs{cs) of (typical) phase points y G Us\Us/2 with C3 < ki{y) < 00 
is dealt with by lemmas 3.36, 3.37, and Corollary 3.38, where it is shown that 
p {Fs{c3)) < C-5, with constants C that can be chosen arbitrarily small by selecting 
the constant C3 1 big enough. The ultimate reason of this measure estimate is 
the following fact: For every point y G Fs{cs) the projection 

n(y) = ^%(^/) e OM 

(where t^^(^y) is the time of the ki{y)-th. collision on the forward orbit of y) will 

have a tubular distance Ztub (^iy)^ ^ Cid from the singularity set S7l~ U STZ'^ , 
where the constant Ci can be made arbitrarily small by choosing the contraction 
coefficients of the powers on the curves p„ arbitrarily small with the 

fei(a) 

help of the result in Appendix II. The upper masure estimate (inside the set dM.) 
of the set of such points U{y) G (Lemma 2 in [S-Ch(1987)]) finally yields the 
required upper bound fx (^5(03)) < C ■ S with arbitrarily small positive constants C 
(if C3 1 is big enough). 
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The listed measure estimates and the obvious fact 

{Us \ Us/2) ^C2-5 

(with some constant C2 > 0, depending only on J) show that there must exist a 
point y & Us\Us/2 with the property (a) above, thus ensuring the sufficiency of the 
point n(y) e J. 

Finally, in the closing section we complete the inductive proof of ergodicity (with 
respect to the number of balls A^) by utilizing Main Lemma 3.5 and earlier results 
from the literature. Actually, a consequence of the Main Lemma will be that excep- 
tional J-manifolds do not exist, and this will imply the fact that no distinct, open 
ergodic components can coexist. 

Appendix I at the end of this paper serves the purpose of making the reading 
of the proof of §3 easier, by providing a chart for the hierarchy of the selection of 
several constants playing a role in the proof of Main Lemma 3.5. 

Appendix II contains a useful (also, potentially useful in the future) uniform 
contraction estimate which is exploited in Section 3. Many ideas of Appendix II 
originate from N. I. Chernov. 



§2. Prerequisites 

Consider the //-dimensional {u > 2), standard, flat torus T*^ = W jTy as the 
vessel containing (> 2) hard balls (spheres) Bi, . . . , Bn with positive masses 
nil, . . . , and (just for simplicity) common radius r > 0. We always assume that 
the radius r > is not too big, so that even the interior of the arising configuration 
space Q (or, equivalently, the phase space) is connected. Denote the center of the 
ball Bi by qi , and let Vi = qi be the velocity of the i-th particle. We investigate 
the uniform motion of the balls Bi, . . . , B^ inside the container with half a unit 

of total kinetic energy: E = - Xli^i — 2' assume that the collisions 

between balls are perfectly elastic. Since — beside the kinetic energy E — the total 
momentum / = rriiVi G is also a trivial first integral of the motion, we 

make the standard reduction / = 0. Due to the apparent translation invariance of 
the arising dynamical system, we factorize the configuration space with respect to 
uniform spatial translations as follows: {qi, . . . , qN) ~ {qi + ci, ■ ■ ■ , qN + ci) for all 
translation vectors a E T'^ . The configuration space Q of the arising flow is then 
the factor torus ^(T^)^ / ~^ ^ 'ji'{N-i) j^in^g i]^q cylinders 

a,j = {(gi, . . . , giv) e T'^(^-i): dist(gi, g,) < 2r} 

(1 < i < J < A^) corresponding to the forbidden overlap between the i-th and j-th 
spheres. Then it is easy to see that the compound configuration point 
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9=(9i,...,9iv)eQ = T-(^-i)\ U Q 



l<i<j<N 



moves in Q uniformly with unit speed and bounces back from the boundaries dCij 
of the cyhnders Cij according to the classical law of geometric optics: the angle of 
reflection equals the angle of incidence. More precisely: the post-collision velocity 
can be obtained from the pre-collision velocity v~ by the orthogonal reflection 
across the tangent hyperplane of the boundary dQ at the point of collision. Here we 
must emphasize that the phrase "orthogonal" should be understood with respect to 
the natural Riemannian metric (the kinetic energy) = X^^^ mi||(igi|p in the 

configuration space Q. For the normalized Liouville measure fu, of the arising flow 
{S^} we obviously have d/j, — const ■ dq ■ dv, where dq is the Riemannian volume 
in Q induced by the above metric, and dv is the surface measure (determined by 
the restriction of the Riemannian metric above) on the unit sphere of compound 
velocities 



The phase space M of the flow {S*} is the unit tangent bundle Q x Ei^~^ of the 
configuration space Q. (We will always use the shorthand notation d = ^{N — 1) 
for the dimension of the billiard table Q.) We must, however, note here that at 
the boundary dQ of Q one has to glue together the pre-collision and post-collision 
velocities in order to form the phase space M, so M is equal to the unit tangent 
bundle Q x S''"^ modulo this identification. 

A bit more detailed definition of hard ball systems with arbitrary masses, as well 
as their role in the family of cylindric billiards, can be found in §4 of [S-Sz(2000)] 
and in §1 of [S-Sz(1999)]. We denote the arising flow by (M, fx). 

In the late 1970s Sinai [Sin(1979)] developed a powerful, three-step strategy for 
proving the (hyperbolic) ergodicity of hard ball systems. This strategy was later 
implemented in a series of papers [K-S-Sz(1989)], [K-S-Sz(1990)-I], [K-S-Sz(1991)], 
and [K-S-Sz(1992)]. First of all, these proofs are inductions on the number of 
balls involved in the problem. Secondly, the induction step itself consists of the 
following three major steps: 

Step I. To prove that every non-singular (i. e. smooth) trajectory segment 
gla,b]^^ with a "combinatorially rich" (in a well deflned sense of Deflnition 3.28 
of [Sim(2002)]) symbolic collision sequence is automatically sufficient (or, in other 
words, "geometrically hyperbolic", see below in this section), provided that the 
phase point xq does not belong to a countable union J of smooth sub-manifolds 
with codimension at least two. (Containing the exceptional phase points.) 
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The exceptional set J featuring this result is negligible in our dynamical consid- 
erations — it is a so called slim set. For the basic properties of slim sets, see again 
below in this section. 

Step II. Assume the induction hypothesis, i. e. that all hard ball systems with 

A^' balls [2 < N' < N) arc (hyperbolic and) ergodic. Prove that there exists a slim 
set i? C M with the following property: For every phase point Xq G M \ i? the 
entire trajectory S^xq contains at most one singularity and its symbolic collision 
sequence is combinatorially rich in the sense of Definition 3.28 of [Sim(2002)], just 
as required by the result of Step I. 

Step III. By using again the induction hypothesis, prove that almost every singular 
trajectory is sufficient in the time interval {to,+oo), where to is the time moment 
of the singular reflection. (Here the phrase "almost every" refers to the volume 
defined by the induced Riemannian metric on the singularity manifolds.) 

We note here that the almost sure sufficiency of the singular trajectories (fea- 
turing Step III) is an essential condition for the proof of the celebrated Theorem 
on Local Ergodicity for semi-dispersive billiards proved by Chernov and Sinai [S- 
Ch(1987)]. Under this assumption, the result of Chernov and Sinai states that in 
any semi-dispersive billiard system a suitable, open neighborhood Uq of any suffi- 
cient phase point G M (with at most one singularity on its trajectory) belongs 
to a single ergodic component of the billiard flow (M, {S'*}teR, //). 

A few years ago Balint, Chernov, Szasz, and Toth [B-Ch-Sz-T(2002)] discovered 
that, in addition, the algebraic nature of the scatterers needs to be assumed, in 
order for the proof of this result to work. Fortunately, systems of hard balls are, by 
nature, automatically algebraic. 

In an inductive proof of ergodicity, steps I and II together ensure that there 
exists an arc-wise connected set C C M with full measure, such that every phase 
point € C is sufficient with at most one singularity on its trajectory. Then the 
cited Theorem on Local Ergodicity (now taking advantage of the result of Step III) 
states that for every phase point G C an open neighborhood Uq of xq belongs 
to one ergodic component of the flow. Finally, the connectedness of the set C and 
/i(M \ C) = imply that the flow (M, {S^}t^s., fj) (now with N balls) is indeed 
ergodic, and actually fully hyperbolic, as well. 

The generator subspace Aij C (1 < ^ < J < A) of the cylinder Cij 

(describing the collisions between the i-th and j-th balls) is given by the equation 

(2.1) Aij = {(gi, ...,qN)e {R'-f : qi = qj], 

see (4.3) in [S-Sz(2000)]. Its ortho-complement Lij C is then deflned by the 
equation 

(2.2) Lij = I (gi, . . . , qn) e {R")^ : qk = ior k ^ i,j, and niiqi + rrijqj = [ , 
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see (4.4) in [S-Sz(2000)]. Easy calculation shows that the cylinder Ci^ (describing 
the overlap of the z-th and j-th balls) is indeed spherical and the radius of its base 

sphere is equal to r^j- = 2r.^^^^^, see §4, especially formula (4.6) in [S-Sz(2000)]. 

The structure lattice C C is clearly the lattice C = (Z^)^ = Z^''. 

Due to the presence of an additional invariant quantity / = X^i^i "^i'^^i; one 
usually makes the reduction J^^^ = and, correspondingly, factorizes the 
configuration space with respect to uniform spatial translations: 

(2.3) (qi, . . . ,qN) ^ iqi + a, . . . ,qN + a), aeT". 

The natural, common tangent space of this reduced configuration space is 



(2.4) Z= Uvi,...,VN)e{Wf : J2m,v, 
[ i=i 

supplied with the inner product 

AT 

{v, v') = ^mi{vi, t;-), 

i=l 

see also (4.1) and (4.2) in [S-Sz(2000)]. 

Collision graphs. Let S^°'-'^^x be a nonsingular, finite trajectory segment with 
the collisions ai, . . . , cx^ listed in time order. (Each ak is an unordered pair (i, j) 
of diff'erent labels i,j G {1, 2, . . . , A^}.) The graph Q = {V,S) with vertex set 
V = {1,2,..., A^} and set of edges S = {ai, . . . , cr^} is called the collision graph 
of the orbit segment S^'^'^^x. For a given positive number C, the collision graph 
G = {V,£) of the orbit segment S^°''''^x will be called C-rich if G contains at least 
C connected, consecutive (i. e. following one after the other in time, according to 
the time-ordering given by the trajectory segment S^'^'^^x) subgraphs. 

Singularities and Trajectory Branches. There are two types of singularities 
of the billiard flow: tangential (or gliding) and double collisions of balls. The first 

means that two balls collide in such a way that their relative velocity vector is 
parallel to their common tangent hyperplane at the point of contact. In this case 
no momentum is exchanged, and on the two sides of this singularity the fiow behaves 
differently: on one side the two balls collide with each other, on the other side they 
fly by each other without interaction. 

The second type of singularity, a double collision means that two pairs of parti- 
cles, and (j, k) {i, j, and k are three different labels) are to collide exactly at 
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the same time. (The case when the two interacting pairs do not have a common 
particle j may be disregarded, since this does not give rise to non-differentiabihty 
of the flow.) 

We are going to briefly describe the discontinuity of the flow {S^} caused by a 
double collisions at time to. Assume first that the pre-coUision velocities of the 
particles are given. What can we say about the possible post-collision velocities? 
Let us perturb the pre-coUision phase point (at time to — 0) infinitesimally, so that 
the collisions at ~ to occur at infinitesimally different moments. By applying the 
collision laws to the arising finite sequence of collisions, we see that the post-collision 
velocities are fully determined by the time-ordered list of the arising collisions. 
Therefore, the collection of all possible time-ordered lists of these collisions gives rise 
to a finite family of continuations of the trajectory beyond to- They are called the 
trajectory branches. It is quite clear that similar statements can be said regarding 
the evolution of a trajectory through a double collision in reverse time. Furthermore, 
it is also obvious that for any given phase point G M there are two, a;-high trees 
and 71 such that (T_) describes all the possible continuations of the positive 
(negative) trajectory S^^'°°^xo {S^~°°'^^xo)- (For the definitions of trees and for 
some of their applications to billiards, cf. the beginning of §5 in [K-S-Sz(1992)].) 
It is also clear that all possible continuations (branches) of the whole trajectory 

Si-oo,oo) 

xo can be uniquely described by all pairs (-B-, -B+) of infinite branches of 
the trees r_ and T+ (S_ c T_, S+ c T+). 

Finally, we note that the trajectory of the phase point a^o has exactly two 
branches, provided that S^xo hits a singularity for a single value t = to, and the 
phase point S^^xq does not lie on the intersection of more than one singularity 
manifolds. In this case we say that the trajectory of xq has a "simple singularity" . 

Other singularities (phase points lying on the intersections of two or more sin- 
gularity manifolds of codimension 1) can be disregarded in our studies, since these 
points lie on a countable collection of codimension- two, smooth submanifolds of the 
phase space, and such a special type of slim set can indeed be safely discarded in 
our proof, see the part Slim sets later in this section. 

Neutral Subspaces, Advance, and Sufficiency. Consider a nonsingular tra- 
jectory segment S^°''^^x. Suppose that a and b are not moments of collision. 

Definition 2.5. The neutral space Mo{S^"-'^^x) of the trajectory segment S^°''''^x at 
time zero (a < < b) is defined by the following formula: 

AfoiS^'^'^h) = {W eZ: 3{S > 0) such that Va e (-6,6) 

V {S"" {Q{x) + aW, V{x))) = ViS^'x) and V {S^ {Q{x) + aW, V{x))) = V{S^x)}. 

[Z is the common tangent space T^Q of the parallelizable manifold Q at any of its 
points q, while V{x) is the velocity component of the phase point x — {Q{x), V{x)).) 
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It is known (see (3) in §3 of [S-Ch (1987)]) that Afo{S^^'^^x) is a hnear subspace 
of Z indeed, and V{x) e J\fo{S^°-'^h). The neutral space J\ft{S^°-'^^x) of the segment 
Sl°-Mx at time t e [a, 6] is defined as follows: 

MtiS^'^'^^x) ^Mo (5["-*'^-*](5*x)) . 

It is clear that the neutral space can be canonically identified with 

A/o(5'["'^lx) by the usual identification of the tangent spaces of Q along the trajectory 
^(-oo,oo)^ (see, for instance, §2 of [K-S-Sz(1990)-I]). 

Our next definition is that of the advance. Consider a non-singular orbit segment 
g[a,b]x with the symbolic collision sequence S = (ui, . . . , (j„), meaning that S^"-'''^x 
has exactly n collisions with dQ, and the z-th collision {1 < i < n) takes place at the 
boundary of the cyhnder Ccr.. For x = (Q, F) G M and W E Z, \\W\\ sufficiently 
smaU, denote Tw{Q, V) := {Q + W,V). 

Definition 2.6. For any 1 <k <n and t e [a,b], the advance 

afc = «((7fe):M(-5l"''lx)^M 

of the collision cr^ is the unique linear extension of the linear functional at = 
a{ak) defined in a sufficiently small neighborhood of the origin of Nt{S^°''^^x) in the 
following way: 

a{ak){W) ■.= tk{x)-tk{S~'TwS'x). 

Here = tk{x) is the time of the fc-th collision ak on the trajectory of x after 
time t = a. The above formula and the notion of the advance functional 

ak = a{ak) : Aft (^t^'^'^x) ^ R 

has two important features: 

(i) If the spatial translation (Q, V) ^ {Q + W, V) {W G Aft (^t"'^Ix)) is carried 
out at time t, then tk changes linearly in and it takes place just ak{W) units of 
time earlier. (This is why it is called "advance".) 

(ii) If the considered reference time t is somewhere between tk-i and t]^, then 
the neutrality of W with respect to ak precisely means that 

W-ak(W)-Vix)eA^„ 

i. e. a neutral (with respect to the collision ak) spatial translation W with the 
advance ak{W) = means that the vector W belongs to the generator space A^-^. 
of the cylinder Ca^. ■ 

It is now time to bring up the basic notion of sufficiency (or, sometimes it is 
also called geometric hyperbolicity) of a trajectory (segment). This is the utmost 
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important necessary condition for the proof of the Theorem on Local Ergodicity for 
semi-dispersive biUiards, [S-Ch(1987)]. 

Definition 2.7. 

(1) The nonsingular trajectory segment S^'^^'^^x (a and b are supposed not to be 
moments of collision) is said to be sufficient if and only if the dimension of 
A/;(5['*'*la;) (t e [a,b]) is minimal, i.e. dim Nt{S^'''^'^x) = 1. 

(2) The trajectory segment S^^-'^^x containing exactly one singularity (a so called 
"simple singularity" , see above) is said to be sufficient if and only if both 
branches of this trajectory segment are sufficient. 

Definition 2.8. The phase point x e M with at most one (simple) singularity is 
said to be sufficient if and only if its whole trajectory S^~°°'°°'^x is sufficient, which 
means, by definition, that some of its bounded segments S^°''''^x are sufficient. 

Note. In this paper the phrase "trajectory (segment) with at most one singularity" 
always means that the sole singularity of the trajectory (segment) , if exists, is simple. 

In the case of an orbit S^~°^'°°^x with at most one singularity, sufficiency means 
that both branches of S^~°^'°°^x are sufficient. 

No accumulation (of collisions) in finite time. By the results of Vaserstein 
[V(1979)], Galperin [G(1981)] and Burago-Ferleger-Kononenko [B-F-K(1998)], in 
any semi-dispersive billiard flow there can only be finitely many collisions in finite 
time intervals, see Theorem 1 in [B-F-K(1998)]. Thus, the dynamics is well defined 
as long as the trajectory does not hit more than one boundary components at the 
same time. 

Slim sets. We are going to summarize the basic properties of codimension-two 
subsets A of a connected, smooth manifold M with a possible boundary and corners. 
Since these subsets A are just those negligible in our dynamical discussions, we shall 
call them slim. As to a broader exposition of the issues, see [E(1978)] or §2 of [K- 
S-Sz(1991)]. 

Note that the dimension dim A of a separable metric space A is one of the three 
classical notions of topological dimension: the covering (Cech-Lebesgue), the small 
inductive (Menger-Urysohn), or the large inductive (Brouwer-Cech) dimension. As 
it is known from general topology, all of them are the same for separable metric 
spaces, see [E(1978)]. 

Definition 2.9. A subset A of M is called slim if and only if A can be covered by 
a countable family of codimension-two (i. e. at least two) closed sets of //-measure 
zero, where /j, is any smooth measure on M. (Cf. Definition 2.12 of [K-S-Sz(1991)].) 
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Property 2.10. The collection of all slim subsets of M is a u-ideal, that is, count- 
able unions of slim sets and arbitrary subsets of slim sets are also slim. 

Proposition 2.11. (Locality). A subset A <Z M is slim if and only if for every 
X E A there exists an open neighborhood U of a; in M such that U (lAis slim. (Cf. 
Lemma 2.14 of [K-S-Sz(1991)].) 

Property 2.12. A closed subset ^ C M is slim if and only if = and 

dim A < dimM - 2. 

Property 2.13. (Integrability) . If A C Mi x M2 is a closed subset of the product 
of two smooth, connected manifolds with possible boundaries and corners, and for 
every x e Mi the set 

Aa, = {ye M2: {x,y) e A} 

is slim in M2, then A is slim in Mi x M2. 

The following propositions characterize the codimension-one and codimension- 
two sets. 

Proposition 2.14. For any closed subset S G M the following three conditions 
are equivalent: 

(i) dimS' < dimM - 2; 

(ii) intS' = and for every open connected set G C M the difference set G\S 
is also connected; 

(iii) intiS" = and for every point a; G M and for any open neighborhood V of x 
in M there exists a smaller open neighborhood W G V of the point x such 
that for every pair of points y, z G W \ S there is a continuous curve 7 in 
the set V\S connecting the points y and z. 

(See Theorem 1.8.13 and Problem 1.8.E of [E(1978)].) 

Proposition 2.15. For any subset S <Z M the condition dim 5" < dimM — 1 is 
equivalent to int^ = 0. (See Theorem 1.8.10 of [E(1978)].) 

We recall an elementary, but important lemma (Lemma 4.15 of [K-S-Sz(1991)]). 
Let A2 be the set of phase points a; G M \ 9M such that the trajectory S^~°^'°°^x 
has more than one singularities (or, its only singularity is not simple). 

Proposition 2.16. The set A2 is a countable union of codimension-two smooth 
submanifolds of M and, being such, is slim. 

The next lemma establishes the most important property of slim sets which 
gives us the fundamental geometric tool to connect the open ergodic components 
of billiard flows. 
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Proposition 2.17. If M is connected, then the complement M\A of a shm F^. set 
^ C M is an arc- wise connected (Gg) set of full measure. (See Property 3 of §4.1 in 
[K-S-Sz(1989)]. The Fa- sets are, by definition, the countable unions of closed sets, 
while the Gs sets are the countable intersections of open sets.) 

The subsets MP and M*. Denote by the set of all phase points a; e M for 
which the trajectory of x encounters infinitely many non-tangential collisions in both 
time directions. The trajectories of the points a; G M \ are lines: the motion is 
linear and uniform, see the appendix of [Sz(1994)]. It is proven in lemmas A. 2.1 and 
A. 2. 2 of [Sz(1994)] that the closed set M \ is a finite union of hyperplanes. It 
is also proven in [Sz(1994)] that, locally, the two sides of a hyper-planar component 
of M \ M'^ can be connected by a positively measured beam of trajectories, hence, 
from the point of view of ergodicity, in this paper it is enough to show that the 
connected components of entirely belong to one ergodic component. This is 
what we are going to do in this paper. 

Denote by M° the set of all phase points x e M''^ the trajectory of which does not 
hit any singularity, and use the notation M-*^ for the set of all phase points x e M# 
whose orbit contains exactly one, simple singularity. According to Proposition 
2.16, the set M# \ (M" U M^) is a countable union of smooth, co dimension-two 
(> 2) submanifolds of M, and, therefore, this set may be discarded in our study of 
ergodicity, please see also the properties of slim sets above. Thus, we will restrict 
our attention to the phase points x e M° U M^. 

The "Chernov-Sinai Ansatz". An essential precondition for the Theorem on 
Local Ergodicity by Chernov and Sinai [S-Ch(1987)] is the so called "Chernov- 
Sinai Ansatz" which we are going to formulate below. Denote by STZ'^ C dM. the 
set of all phase points xq = (qo^vo) G corresponding to singular reflections 
(a tangential or a double collision at time zero) supplied with the post-collision 
(outgoing) velocity vq. It is well known that STZ'^ is a compact cell complex with 
dimension 2d — 3 = dimM — 2. It is also known (see Lemma 4.1 in [K-S-Sz(1990)-I], 
in conjunction with Proposition 2.16 above) that for i^i-almost every phase point 
Xq G STZ'^ the forward orbit S^'^'^^^xq docs not hit any further singularity. (Here 
ui is the Riemannian volume of STZ'^ induced by the restriction of the natural 
Riemannian metric of M.) The Chernov-Sinai Ansatz postulates that for z/i-almost 
every xq G Sn+ the forward orbit S^°'°°^ Xq is sufficient (geometrically hyperbolic). 

The Theorem on Local Ergodicity. The Theorem on Local Ergodicity for 
semi-dispersive billiards (Theorem 5 of [S-Ch(1987)]) claims the following: Let 
(M, {S'*}t£R, /x) be a semi-dispersive billiard flow with the property that the smooth 
components of the boundary dQ of the configuration space are algebraic hyper- 
surfaces. (The cylindric billiards automatically fulfill this algebraicity condition.) 
Assume - further - that the Chernov-Sinai Ansatz holds true, and a phase point 
xo G (M° U M^) \ is sufficient. 
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Then some open neighborhood t/o C M of xq belongs to a single ergodic compo- 
nent of the flow (M, {S'*}teR, /i). (Modulo the zero sets, of course.) 

§3. The Exceptional J-Manifolds 
(The asymptotic measure estimates) 



First of all, we define the fundamental object for the proof of our theorem. 

Definition 3.1. A smooth submanifold J C intM of the interior of the phase space 
M is called an exceptional J -manifold (or simply an exceptional manifold) with a 
negative Lyapunov function Q if 

(1) dim J = 2d-2{= dimM - 1); 

(2) the pair of manifolds (J, dJ) is diffeomorphic to the standard pair 

where B^^~'^ is the closed unit ball of E^'^"^; 

(3) J is locally flow-invariant, i. e. Vx G J 3a{x), b{x), a{x) < < b{x), such 
that S^x e J for aU t with a{x) < t < b{x), and 5"(^)a; e dJ, S^^^'^x e dJ; 

(4) the manifold J has some thin, open, tubular neighborhood Uq in intM, and 
there exists a number T > such that 

(i) S'^ ^ ^"'^ ^ ^' orbit segments {x e C/q) are non-singular, 
hence they share the same symbolic collision sequence E; 

(ii) Vx G Uq the orbit segment S^^^'^^x is sufficient if and only if x ^ J; 

(5) \/x E J we have Q{n{x)) := {z{x), w{x)) < —ci < for a unit normal vector 
fleld n{x) = {z{x), w{x)) of J with a fixed constant ci > 0; 

(6) the set W of phase points a; e J never again returning to J (After first leaving 
it, of course. Keep in mind that J is locally flow-invariant) has relative measure 

greater than 1 — 10~^ in J, i. e. . ^. > 1 — 10~^, where ui is the hypersurface 

measure of the smooth manifold J. 

Remark. The above definition is, by nature, fairly technical, thus a short com- 
menting of it is due here. Once we make the induction hypothesis, i. e. we assume 
that the (hyperbolic) ergodicity and the Chernov-Sinai Ansatz hold true for any 
hard ball systems with less than N balls (regardless of the masses and the ra- 
dius r), the only way for two distinct ergodic components to co-exist is when they 
arc separated by an exceptional manifold J described in the above definition. This 
is proved in §4 below. 
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We begin with an important proposition on the structure of forward orbits 
^[o,oo)^ for xeJ. 

Proposition 3.2. For //i-almost every x E J the forward orbit S^^'°°''x is non- 
singular. 

Proof. According to Proposition 7.12 of [Sini(2003)], the set 



Jn 



.t>0 



of forward singular points x & J is a countable union of smooth, proper submanifolds 
of J, hence it has //i-measure zero. □ 

In the future we will need 

Lemma 3.3. The concave, local orthogonal manifolds passing through points 
y & J are uniformly transversal to J. 

Note. A local orthogonal manifold S C intM is obtained from a codimension-one, 
smooth submanifold Ei of intQ by supplying Si with a selected field of unit normal 
vectors as velocities. E is said to be concave if the second fundamental form of Ei 
(with respect to the selected field of normal vectors) is negative semi-definite at 
every point of Ei. Similarly, the convexity of E requires positive semi-definiteness 
here, see also §2 of [K-S-Sz(1990)-I]. 

Proof. We will only prove the transversality. It will be clear from the uniformity 
of the estimates used in the proof that the claimed transversalities are actually 
uniform across J. 

Assume, to the contrary of the transversality, that a concave, local orthogonal 
manifold E(|/) is tangent to J at some y E J. Let {Sq, B5q) be any vector of TyNl 
tangent to E(|/) at y. Here S < is the second fundamental form of the projection 
qiT^{y)) = Ei(y) of E(y) at the point q = q{y). The assumed tangency means 
that {5q, z) + {B5q^ w) = 0, where n{y) = {z{y), w{y)) = {z,w) is the unit normal 
vector of J at y. We get that {dq, z + Bw) = for any vector Sq G v{y)-^. We note 
that the components z and w oi n are necessarily orthogonal to the velocity v{y), 
because the manifold J is locally flow-invariant and the velocity is normalized to 1 
in the phase space M. The last equation means that z = —Bw, thus Q{n{y)) = 
{z, w) = {—Bw, w) > 0, contradicting to the assumption Q{n{y)) < — ci of (5) in 
3.1. This finishes the proof of the lemma. □ 

In order to formulate the main result of this section, we need to define two 
important subsets of J. 

Definition 3.4. Let 

A = |x e J| S^^'^^^x is nonsingular and dim^Vo (^S^^'°°^x^ ^ ^} ' 
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B=i^xeJ\ S^°'°°^x is nonsingular and dimjVo (^5[°'°°)a;j > l| . 

The two Borel subsets A and B of J are disjoint and, according to Proposition 
3.2 above, their union AU B has fuU /xi-measure in J. 
The anticipated main result of this section is 

Main Lemma 3.5. Use all of the above definitions and notations. We claim that 
A^$. 

Proof. The proof will be a proof by contradiction, and it will be subdivided into 
several lemmas. Thus, from now on, we assume that A = $. 

First, select and fix a non-periodic point (a "base point") xq G B. Following the 
fundamental construction of local stable invariant manifolds [S-Ch(1987)] (see also 
§5 of [K-S-Sz(1990)-I]), for any j/ e M and any t > we define the concave, local 
orthogonal manifolds 



Eliy) = SCy, ({(g, v{yt)) e M| g - ± v{yt)} \ {S, U , 

^liy) = scy [s-'^liv)] , 

where 5i := G M| Tx G STZ~^ (the set of phase points on singularities of order 
1), S-i :— {xgMI — X e Si} (the set of phase points on singularities of order 
— 1), yt = S*y, and SCy{ . ) stands for taking the smooth component of the given 
set that contains the point y. The local, stable invariant manifold 7*(y) of y is 
known to be a superset of the C^-limiting manifold lim^^oo '^oiv)- 

For any j/ G M we use the traditional notations 



T{y) = min {t > 0| S^y G ^M} , 
^^'^^ T{y) = S^y^y 

for the first hitting of the collision space 9M. The first return map (Poincare 
section, collision map) T : 9M 9M (the restriction of the above T to dNl) is 
known to preserve the finite measure u that can be obtained from the Liouville 
measure by projecting the latter one onto dM. along the flow. Following 4. of 
[K-S-Sz(1990)-II], for any point y G intM (with T{y) < oo, T{—y) < oo, where 
~y = {Qj ~'>j) for y — {Qj '>j)) we denote by ztub{y) the supremum of all radii p > 
of tubular neighborhoods Vp of the projected segment 

q{{S'y\ -T{-y)<t<T{y)])^Cl 
for which even the closure of the set 
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{{q, viy))eM\qeV,} 

does not intersect the set Si U S-i. We remind the reader that both Lemma 2 of 
[S-Ch(1987)] and Lemma 4.10 of [K-S-Sz(1990)-I] use this tubular distance function 
ztubi-) (despite the notation z{.) in those papers), see the important note 4. in 
[K-S-Sz(1990)-II]. 

On all the arising local orthogonal manifolds, appearing in the proof, we will 
always use the so called 5g-metric to measure distances. The length of a smooth 
curve with respect to this metric is the integral of \\5q\ \ along the curve. The proof 
of the Theorem on Local Ergodicity [S-Ch(1987)] shows that the Sq-metric is the 
relevant notion of distance on the local orthogonal manifolds S, also being in good 
harmony with the tubular distance function ztub{ ■ ) defined above. 

The tangent vector u{xo) = [Sqo, Svq) is defined as follows: 

For a large constant Lq ^ 1 (to be specified later), wc select first a non-collision 
time £3 ^ 1 in the following way: Thanks to our hypothesis (5) in 3.1 and the hy- 
persurface measure expansion theorem of [Ch-Sim(2006)], the hypersurface measure 
of S*{J) grows at least linearly in t, as t — > 00. As a consequence, the distances 
between J and nearby points will shrink at least linearly in t, as t ^ 00. The 
theorem of Appendix II below claims that for a large enough (non-collision) time 
63 3> 1 the phase point — S'^^xq has a unit tangent vector {6qo,dvo) e E^{xc^) 
such that the normalized tangent vector 



(3.8) .(.0) = {5qo, 5vo) := ||^^^-.3)(^,^,^,^)|| ^ ^ (-0) 

is transversal to J, and the expansion estimate 

||r~ II > ^^0 
ll%3ll 

holds true or, equivalently, we have the contraction estimate 



(3 9) Il^^g3ll _±_ 

^ ^ ^ ||<^go|| 2Lo' 

where 

{5qo,5vo) 



\\{DS-^^){6qo,6vo] 
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Remark. Almost every phase point xq of the hypersurface J satisfies the hypothe- 
ses of Appendix II (on the connected coUision graphs). This is indeed so, since 
the proof of Theorem 6.1 of [Sim(1992)-I] works without any essential change not 
only for singular phase points, but also for the points of the considered exceptional 
manifold J. The only important ingredient of that proof is the transversality of 
the spaces E^{x) to J, provided by Lemma 3.3 above. According to that result, 
typical phase points x E J (with respect to the hypersurface measure of J) indeed 
enjoy the above property of having infinitely many consecutive, connected collision 
graphs on their forward orbit S^^'^^^xq. 

We choose the orientation of the unit normal field n{x) {x G J) of J in such a 
way that {n{xo), {Sqo, Svq)) < 0, and define the one-sided tubular neighborhood Us 
of radius 5 > as the set of all phase points 7a; (s), where x E J, < s < 6. Here 
7a;( . ) is the geodesic line passing through x (at time zero) with the initial velocity 
n{x), X E J. The radius (thickness) 5 > here is a variable, which will eventually 
tend to zero. We are interested in getting useful asymptotic estimates for certain 
subsets of Us, as 5 ^ 0. 

Our main working domain will be the set 

Do ^ eU5^\j\y ^ [J S'^ {STZ-) , 3 a sequence 
(3.10) *>o 

tn y oo such that 5*"?/ eUso\J, n = 1, 2, . . . |, 

a set of full |U-measure in Us^- We will use the shorthand notation Uq = Us^ for a 
fixed, small value Sq. 

On any manifold So(2/) we define the smooth field Xy^t{y') {y' £ ^oiv) ^ Uq) 
of unit tangent vectors of S^(y) fl Uq as follows: 



(3.11) x,4y') = ^y^'-' ^^'^^ 



where 11^^^^^/ denotes the orthogonal projection of M*^ © M'^ onto the tangent space 
of Ilo(y) at the point y' e Eo(y) n Uq. 

If, in the construction of the manifolds TiQ{y), the time t is large enough, i. e. 
t > C3 for a suitably large constant C3 ^ £3, the points y, y' are close enough to xq, 
and y' e So(j/), then the tangent space Ty'T,Q{y) will be close enough to the tangent 
space TcQ7*(a;o) of the local stable manifold 7*(xo) of xq, so that the projected copy 
'^y,t{y') of {Sqo, 6vo) (featuring (3.11)) will undergo a contraction by a factor of at 
least Lq ^ between time and 63, let alone between time and C3, that is. 
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for all t > C3. We note that the tangent space %;gY{^o) of the local stable manifold 
7*(xo) makes sense, even if the latter object does not exist: this tangent space can 
be obtained as the positive subspace of the operator B{xq) defined by the continued 
fraction (2) in [S-Ch(1987)] or, equivalently, as the intersection of the inverse images 
of stable cones of remote phase points on the forward orbit of xq. All the necessary 
upper estimates for the mentioned angles between the considered tangent spaces 
follow from the well known result stating that the difference (in norm) between the 
second fundamental forms of the S'*-images {t > 0) of two local, convex orthogonal 
manifolds is at most see, for instance, inequality (4) in [Ch(1982)]. These facts 
imply, in particular, that the vector in the numerator of (3.11) is actually very close 
to (Sqo, Svq), thus its magnitude is almost one. 

For any y E Dq let — tk{y) {0 < ti < t2 < ■ ■ ■) be the time of the k-th collision 
CTfc on the forward orbit S^^'°°^y of y. Assume that the time t in the construction 
of Sq(j/) and Xy^t is between cTfc_i and Ufc, i. e. tk-i{y) < t < tk{y)- We define the 
smooth curve py,t = Py.t{s) (with the arc length parametrization s, < s < h{y,t)) 
as the maximal integral curve of the vector field Xy^t emanating from y and not 
intersecting any forward singularity of order < k, i. e. 



(3.13) 



Py,t(0) = y, 

Py,t{ ■ ) does not intersect any singularity of order < /c, 
Py^t is maximal among all curves with the above properties. 



We remind the reader that a phase point x lies on a singularity of order k {k eN) 
if and only if the A;-th collision on the forward orbit S^^'°°'>x is a singular one. 
It is also worth noting here that, as it immediately follows from definition (3.6) 
and (3.13), the curve py t can only terminate at a boundary point of the manifold 

my)nUo. 

Remark 3.14. From now on, we will use the notations Eq^/), Xy^k, and py^k for 
^o'iy)^ ^y,tl^ and py,^*, respectively, where % = tl{y) = \ {tk-i{y) + tk{y))- 

Due to these circumstances, the curves Py,ti = Py,k can now terminate at a point 
z such that z is not on any singularity of order at most k and S^i^z is a boundary 
point of E^J(y), so that at the point S^i^z the manifold E^5(y) touches the boundary 
of the phase space in a nonsingular way. This means that, when we continuously 
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move the points Py,k{s) by varying the parameter s between and h{y, k), either 
the time tk iPy,k{s)) or the time tk-i {Py,k{s)) becomes equal to t^. = t^iy) when the 
parameter vahic s reaches its maximum value h(y, k). The length of the curve py^k is 
at most Sq, and an elementary geometric argument shows that the time of collision 
tk{Py,k{s)) (or tk-i{py,k{s))) can only change by at most the amount of c*v^, as 
s varies between and h{y,k). (Here c* is an absolute constant.) Thus, we get 
that the unpleasant situation mentioned above can only occur when the difference 
tk{y) — tk-i{y) is at most c*y/5o. These collisions have to be and will be excluded as 
stopping times k2{y), f-iiv) and ki{y) in the proof below. Still, everything works by 
the main result of [B-F-K(1998)], which states that there is a large positive integer 
no and a small number f3 > such that amongst any collection of no consecutive 
collisions there are always two neighboring ones separated from each other (in time) 
by at least /?. Taking c* < P shows that the badly behaved collisions - described 
above - can indeed be excluded from our construction. 

As far as the terminal point py^k{h{y, k)) of py^k is concerned, there are exactly 
three, mutually exclusive possiblities for this point: 

(A) py^k{h{y,k)) G J and this terminal point does not belong to any forward 
singularity of order < /c, 

(B) py^k{h{y, k)) lies on a forward singularity of order < k, 

(C) the terminal point py^k{h{y, k)) does not lie on any singularity of order < k 
but lies on the part of the boundary OUq of Uq different from J = J x {0} and 
J X {5o}. 

Remark 3.15. Under the canonical identification Uq = J x [0, 5o) of C/q via the 
geodesic lines perpendicular to </, the above mentioned part of dUo (the "side" of 
Uq) corresponds to dJ x [0, Sq). Therefore, the set of points with property (C) 
inside a layer Us {S < Sq) will have //-measure o{S) (actually, of order 5^), and 
this set will be negligible in our asymptotic measure estimations, as 5 ^ 0. The 
reason why these sets are negligible, is that in the indirect proof of Main Lemma 
3.5, a contradiction will be obtained (at the end of §3) by comparing the measures 
of certain sets, whose measures are of order const • 5. That is why in the future we 
will not be dealing with any phase point with property (C). 

Should (B) occur for some value of A; (A; > 2), the minimum of all such integers 
k will be denoted hj k = k{y). The exact order of the forward singularity on which 
the terminal point p^^ {h{y, k)) lies is denoted by ki = ki{y) (< k{y)). If (B) does 

not occur for any value of k, then we take k{y) = ki{y) = oo. 

We can assume that the manifold J and its one-sided tubular neighborhood 
Uq = Uso are already so small that for any y e Uq no singularity of S^^'°°^y can 
take place at the first collision, so the indices k and ki above are automatically at 
least 2. For our purposes the important index will be ki = ki{y) for phase points 
y e ^0- 
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Remark 3.16. Refinement of the construction. Instead of selecting a single 
contracting unit vector {5qo, Svq) in (3.8), we should do the following: Choose a 
compact set Kq C B with the property 

^i(j) 

Now the running point x € Kq will play the role of xq in the construction of the 
contracting unit tangent vector u{x) := {dqo, Svq) G E^{x) on the left-hand-side of 
(3.8). For every x & Kq there is a small, open ball neighborhood B{x) of x and a big 
threshold c^{x) ^ 1 such that the contraction estimate (3.9) and the transversality 
to J hold true for u{y) and £3 = Cs{x) for all y G B{x). 

Exactly the same way as earlier, one can also achieve that the weaker contraction 
estimate Lq^ of (3.12) holds true not only for t > and u{x), but also for any 
projected copy of it appearing in (3.11) and (3.12), provided that y, y' G B{x), and 
t > C3{x). 

Now select a finite subcover \Jhi=i B{xi) of Kq, and replace J by Ji = J fl 
[j^^i B{xi), Us by Ug = Us ^ [j^^i B{xi) (for d < Sq) and, finally, choose the 
threshold C3 to be the maximum of all thresholds cs{xi) for z = 1, 2, . . . , n. In this 
way the assertion of Corollary 3.18 below will be true. 

We note that the new exceptional manifold Ji is no longer so nicely "round 
shaped" as J, but it is still pretty well shaped, being a domain in J with a piecewise 
smooth boundary. 

The reason why we cannot switch completely to a round and much smaller man- 
ifold B{x) f] J is that the measure //i(J) should be kept bounded from below after 
having fixed Lq, see the requirement 3 in Appendix I. 

Remark 3.17. When defining the returns of a forward orbit to Us, we used to 
say that "before every new return the orbit must first leave the set C/5". Since the 
newly obtained J is no longer round shaped as it used to be, the above phrase is 
not satisfactory any longer. Instead, one should say that the orbit leaves even the 
^^-neighborhood of Us, where k is two times the diameter of the original J. This 
guarantees that not only the new Us, but also the original Us will be left by the 
orbit, so we indeed are dealing with a genuine return. This note also applies to two 
more shrinkings of J that will take place later in the proof. 

In addition, it should be noted that, when constructing the vector field in (3.11) 
and the curves py^t, an appropriate directing vector u{xi) needs to be chosen for 
(3.11). To be definite and not arbitrary, a convenient choice is the first index 
i G {1,2, . . .,n} for which y G B{xi). In that way the whole curve py^t will stay 
in the slightly enlarged ball B'{xi) with double the radius of B{xi), and one can 
organize things in such a way that the required contraction estimates of (3.12) be 
still true, even in these enlarged balls. 

In the future, a bit sloppily, Ji will be denoted by J, and U'^ by Us- 
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As an immediate corollary of (3.12), the uniform transvarsality of the field 
Xy^t{y') to J and Remark 3.16, we get 

Corollary 3.18. For the given sets J, Uq, and the large constant Lq we can select 
the threshold C3 > large enough so that for any point y E Dq any time t with 
C3 < t < %^(j/)(j/) the (^Qf-expansion rate of between the curves Py^y^ and 

(Py;k{y)) ^^^^ -^0 i- ^- for any tangent vector {6qo, 6vo) of Pj^ fc(y) we 
have 

\\Sqt 



\\Sqo\ 

where {5qt, 5vt) = {DS^){5qo, Svq). 

Remark 3.19. The reason why there is no expansion from time C3 until time t is 
that all the image curves S'^ (^Pyk{y)) (^3 < "T" < t) are concave, according to the 
construction of the curve p -rt \- 

An immediate consequence of the previous result is 

Corollary 3.20. For any y G Dq with k{y) < 00 and t-j^_^^y^_-^{y) > C3, and for any 
t with t^^^y^_^{y) < t < %^y^{y), we have 

(3.21) ztub (S'y) < {PyMy)) < •^)' 

where lq (^p^ ^y-^^ denotes the ^qf-length of the curve p^ and C4 > is a constant, 

independent of Lq or C3, depending only on the (asymptotic) angles between the 
curves p -r, s and J. 

Proof. The manifold J and the curves p^ are uniformly transversal, as it follows 

immediately from the uniformity of the transversality of the field Xy^tiv') to J. This 
is why the above constant C4, independently of Lq, exists. □ 

By further shrinking the exceptional manifold J a little bit, and by selecting a 
suitably thin, one-sided neighborhood Ui = Us^ of J, we can achieve that the open 
25i-neighborhood of Ui (on the same side of J as Uq and Ui) is a subset of Uq. 

For a varying 5, < 5 < 5i, we introduce the layer 



Us = \y ^ {UsX Us/2) riDo\3 a sequence tn 00 
(3.22) ^ 

such that G {Us \ Us/2) for all n L 
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Since almost every point of the layer (t/^ \ H-Dq returns infinitely often to this 
set and the asymptotic equation 



(3.23) {{Us \ Us/2) n Do) - -//I (J) 



holds true, we get the asymptotic equation 



(3.24) //(C/,)~^/xi(J). 
We will need the following subsets oi Us'- 



Us{oo) = {2/ e Us\ ki{y) = 00} . 

Here C3 is the constant from Corollary 3.18, the exact value of which will be specified 
later, at the end of the proof of Main Lemma 3.5. Note that in the first line of (3.25) 
the case ki{y) = cxo is included. By selecting the pair of sets (C/i, J) small enough, 
we can assume that 



(3.26) ztub{y) > C4(^i Vy e t/i. 

This inequality guarantees that the collision time 'tki(y)^y^ (^ ^ cannot be near 
any return time of y to the layer {Us \ Us/2), for 6 < di, provided that y G Us{c3). 
More precisely, the whole orbit segment S^~'^^~^^''^^^'>^z will be disjoint from Ui, 
where z = S'y, %^(^)_i(y) < t < H^^y^)- 

Let us consider now the points y of the set Us{oo). We observe that for any point 
y e Us{oo) the curves Py,k{s) (0 < s < h{y,k)) have a C^-limiting curve Py, oo{s) 
(0 < s < h{y, 00)), with h{y, k) h{y, 00), as A; — > 00. 

Indeed, besides the concave, local orthogonal manifolds Sq(j/) = Sq'=(j/) of (3.6) 
(where tl = t\{y) = \ {tk-i{y) + tk{y)))-, let us also consider another type of con- 
cave, local orthogonal manifolds defined by the formula 

tl{y) = E*^(2/) = SCy (5-*^ (^SCy^. {y' e M\ q{y') = g(2/t*)})) , 

the so called "candle manifolds", containing the phase point y G Us{oo) in their 
interior. It was proved in §3 of [Ch(1982)] that the second fundamental forms 
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B (Eq(|/), y) < are monotone non-increasing in fc, while the second fundamental 
forms B (^^{y)^ < are monotone increasing in /c, so that 

B(tl{y),y)<B{^l{y),y) 
is always true. It is also proved in §3 of [Ch(1982)] that 

lim B (Sg(y), y) = lim B y) := B^ < 

uniformly in y, and these two-sided, monotone curvature limits give rise to uniform 
C^-convergences 

lim SS(y) = hm = S-(y), 

I— >oo t— >oo 

and the limiting manifold Eo°(2/) ^^e local stable invariant manifold 7*(j/) of y, 
once it contains y in its smooth part. These monotone, two-sided limit relations, 
together with the definition of the curves Py,ti = Py,k prove the existence of the 
C^-limiting curve Py,oo = hnifc^oo Pi/,A;, h{y,k) — > h{y, oo), as k ^ oo. They also 
prove the inclusion Py, oo ([0, h{y, oo)]) C 7^(y). 

Lemma 3.27. p, (Us \ Usic^)) = o(5), as 5 ^ 0. 

Proof. The points y of the set Us \ Usics) have the property t-^^(^y-^_i{y) < C3. 
By doing another slight shrinking to J, the same way as in Remark 3.16, we can 
achieve that %^(j,)(l/) < 2c3 for all y G Us\U s{c^), < 5 < 5i. This means that the 
terminal point n(|/) = Pyki{y) {^{Vi ^)) curve Py'j^-^(^y^ lies on the singularity 

set 

U s-'{sn-), 

0<t<2c3 

hence all points of the set Us \ Us{cs) are at most at the distance of S from the 
singularity set mentioned above. 

This singularity set is a compact collection of co dimension-one, smooth subman- 
ifolds (with boundaries), each of which is uniformly transversal to the manifold J. 
This uniform transversality follows from Lemma 3.3 above, and from the fact that 
the inverse images S~*{STl~) (t > 0) of singularities can be smoothly foliated with 
local, concave orthogonal manifolds. Thus, the (5-neighborhood of this singularity 
set inside Us clearly has |U-measure o(5), actually, of order < const ■6'^. □ 

For any point y G Us{oo) we define the return time ^2 = ^2(y) as the infimum of 
all numbers ^2 > C3 for which there exists another number ti, < ti < ^2, such that 
S^^y ^ Uq and S*^{y) G (Us \ Us/2) H Dq. Let k2 = k2{y) be the unique natural 
number for which tk^-iiy) < hiy) < tk^^iy)- 
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Lemma 3.28. For any point y e Us^oo) the projection 

n(y) := PvMiv) (Hy^My))) 

is a forward singular point of J. 

Proof. Assume that the forward orbit of Ii{y) is non-singular. The distance 
dist(/S'*2j/, J) between S^^y and J is bigger than 5/2. According to the contrac- 
tion result 3.18, if the contraction factor Lq^ is chosen small enough, the distance 

between 5"*^ (n(y)) and J stays bigger than 5/4, so S^^ (n(y)) E Uo\J will be true. 
This means, on the other hand, that the forward orbit of n(y) is sufficient, accord- 
ing to (4)/(ii) of Definition 3.1. However, this is impossible, due to our standing 
assumption A = $. □ 

Lemma 3.29. The set Us{oo) is actually empty. 

Proof. Just observe that in the previous proof the whole curve Py,k2{y) can be 
slightly perturbed (in the C°° topology, for example), so that the perturbed curve 
Py emanates from y and terminates on a non-singular point n(y) of J (near Il{y)), so 
that the curve py still "lifts" the point n(j/) up to the set (Ug \ Us 74) flDo if we apply 
5*2. This proves the existence of a non-singular, sufficient phase point Ii{y) e A, 
which is impossible by our standing assumption A = 0. Hence Us{oo) = 0. □ 

Next we need a useful upper estimate for the /i-measure of the set Usics) as 
5—^0. We will classify the points y E Us{cs) according to whether S^y returns to 
the layer (Us \ Us/2) H-Dq (after first leaving it, of course) before the time t-^^^y^_^{ii) 
or not. Thus, we define the sets 

Esicz) ^{y e Us{cs) \ 3 < ti < t2 < h^^y)_M 
(3.30) such that S'^y ^Uo, S'^ye {Us \ Us/2) n Do}, 

Fsics) =Us{c3)\Esic3). 

Recall that the threshold t-^^^y^_^{y), being a collision time, is far from any possible 

return time t2 to the layer {Us \ Us/2) H Dq, see the remark right after (3.26). 

Now we will be doing the "slight shrinking" trick of Remark 3.16 the third (and 
last) time. We slightly further decrease J to obtain a smaller Ji with almost the 
same /Ui-measure. Indeed, by using property (6) of 3.1, inside the set J H B we 
choose a compact set Ki for which 

> 1 - 10-6 

and no point of Ki ever returns to J. For each point x G Ki the distance between 
the orbit segment S^"'°''^^^x and J is at least e{x) > 0. Here is needed to guarantee 
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that we certainly drop the initial part of the orbit, which still stays near J, and 
C3 was chosen earlier. By the non-singularity of the orbit segment S^^°''^^^x and by 
continuity, the point x G Ki has an open ball neighborhood B{x) of radius r{x) > 
such that for every y G B{x) the orbit segment S^'^"''^^^y is non-singular and stays 
away from J by at least e{x)/2. Choose a finite covering UILi ^i^i) ^ -^1 -^i' 
replace J and Us hj their intersections with the above union (the same way as it 
was done in Remark 3.16), and fix the threshold value of 5i so that 

5i < ^ min{e(a;i)| z = 1, 2, . . . , n}. 

In the future we again keep the old notations J and Us for these intersections. In 
this way we achieve that the following statement be true: 

any return time ^2 of any point y G (Us \ Us/2) H Dq to 
(Us \ Us/2) n -Do is always greater than ior < 5 < 5i. 

Just as in the paragraph before Lemma 3.28, for any phase point y G Es{c3) we 
define the return time ^2 = ^2(2/) as the infimum of all the return times t2 of y 
featuring (3.30). By using this definition of hiy), formulas (3.30)-(3.31), and the 
contraction result 3.18, we easily get 

Lemma 3.32. If the contraction coefficient Lq^ in 3.18 is chosen suitably small, 
for any point y G Es{cz) the projected point 




(3.33) Il{y) := p^-^^^y^{h{yMy))) ^ J 

is a forward singular point of J. 

Proof. Since ^2(2/) < %^(y)_i(2/)7 we get that, indeed, n(j/) G J. Assume that the 
forward orbit of n(y) is non-singular. 

Since S'^^^y'^y G (Us \ Us/2) n -Dq, we obtain that dist [s'^^^y'^y, > 5/2. On 
the other hand, by using (3.31) and Corollary 3.18, we get that for a small enough 
contraction coefficient Lq ^ the distance between S^^^y^y and S^'^'^y'^ (n(2/)) is less 
than 5/4:. (The argument is the same as in the proof of Lemma 3.28.) In this way 
we obtain that S'*^^^) (U(y)) e Uo\J, so Il(y) G A, according to condition (4)/(u) 
in 3.1, thus contradicting to our standing assumption A = 0. This proves that, 
indeed, 11 (y) is a forward singular point of J. □ 

Lemma 3.34. The set Es(cs) is actually empty. 
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Proof. The proof will be analogous with the proof of Lemma 3.29 above. Indeed, 
we observe that in the previous proof for any point y e Esics) the curve Pyt2{y) 
can be slightly perturbed (in the C°° topology), so that the perturbed curve py 
emanates from y and terminates on a non-singular point n(y) of J, so that the 
curve /5y still "lifts" the point n(y) up to the set (Us \ Us/2) H Dq if we apply S^'^. 
This means, however, that the terminal point n(y) of Py is an element of the set 
A, violating our standing assumption A — fj). This proves that no point y G £'5(03) 
exists. □ 

For the points y G ^5(03) = Usics) we define the projection Il{y) by the formula 

(3.35) n(y) := e dM. 

Now we prove 

Lemma 3.36. For the measure v (l[{Fs{c2))j of the projected set il{Fs{c3)) C 
dM. we have the upper estimate 

iy(u{Fs{c3))') <C2C^L^'d, 

where C2 > is the geometric constant (also denoted by C2) in Lemma 2 of [S- 
Ch(1987)] or in Lemma 4.10 of [K-S-Sz(1990)-I], C4 is the constant in (3.21) above, 
and u is the natural T-invariant measure on that can be obtained by projecting 
the Liouville measure p onto dM. along the billiard flow. 

Proof. Let y G ^5(03). From the inequality t^_^(^y^_i{y) > C3 and from Corollary 

3.20 we conclude that Zfub (j^iy)j < c^Lq^S. This inequality, along with the fun- 
damental measure estimate of Lemma 2 of [S-Ch(1987)] (see also Lemma 4.10 in 
[K-S-Sz(1990)-I]) yield the required upper estimate for u (fl{Fs{c3))^ . □ 

The next lemma claims that the projection 11: ^5(03) — > dM (considered here 
only on the set Fs{cs) = Usic^)) is "essentially one-to-one", from the point of view 
of the Poincare section. 

Lemma 3.37. Suppose that yi, 2/2 G Fs{cs) are non-periodic points (S < Si), and 
n(yi) = n(y2)- We claim that yi and y2 belong to an orbit segment S of the 
billiard flow lying entirely in the one-sided neighborhood Uq of J and, consequently, 
the length of the segment S is at most l.ldiam( J). 

Remark. We note that in the length estimate l.ldiam(J) above, the coefficient 
1.1 could be replaced by any number bigger than 1, provided that the parameter 
5 > is small enough. 
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Proof. The relation Il{yi) = ^{1/2) implies that yi and y2 belong to the same 
orbit, so we can assume, for example, that y2 = S°'yi with some a > 0. We need 
to prove that S'l'^'^'yi C Uq- Assume the opposite, i. e. that there is a number ti, 
< ti < a, such that S^^yi ^ Uq. This, and the relation S'^yi E (Us \ Us/2) H Dq 
mean that the first return of yi to (Us \ Us/2) ^ Do occurs not later than at time 
t = a. On the other hand, since n(|/i) = n(5'"2/i) and yi is non-periodic, we get that 
> a, see (3.35). The obtained inequality ^^(^^^.^(yi) > a > ^2(2/1), 
however, contradicts to the definition of the set ^5(03), to which yi belongs as 
an element, see (3.30). The upper estimate l.ldiam(J) for the length of S is an 
immediate corollary of the containment S C Uq. □ 

As a direct consequence of lemmas 3.36 and 3.37, we obtain 

Corollary 3.38. For all small enough 5 > 0, the inequality 

At(-^<5(c3)) < I.IC2C4LQ ^Miam(J) 

holds true. 

Finishing the Indirect Proof of Main Lemma 3.5. 

It follows immediately from Lemma 3.27 and corollaries 3.29, 3.34, and 3.38 that 

l^iPs) < 1.2c2C4diam(J)LQ ""^(J 

for all small enough 6 > 0. This fact, however, contradicts to (3.24) if Lq^ is 
selected so small that 

1.2c2C4diam(J)Lg ^ < ^/ii(J*), 

where J* stands for the original exceptional manifold before the three slight shrink- 
ings in the style of Remark 3.16. Clearly, /ii(J) > (1 - 10~^)/ii(J*). The obtained 
contradiction finishes the indirect proof of Main Lemma 3.5. □ 

§4. Proof of Ergodicity 
Induction on N 

By using several results of Sinai [Sin(1970)], Chernov-Sinai [S-Ch(1987)], and 
Kramli-Simanyi-Szasz, in this section we finally prove the ergodicity (hence also 
the Bernoulli property; see Sinai's results in [Sin(1968)], [Sin(1970)], and [Sin(1979)] 
for the K-mixing property, and then by Chernov-Haskell [C-H(1996)] and Ornstein- 
Weiss [0-W(1998)] the Bernoulli property follows from mixing) for every hard ball 
system (M, {5**}, //), by carrying out an induction on the number N (> 2) of 
interacting balls. 
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In order for the proof to work, from now on we assume that the Chernov-Sinai 
Ansatz is true for every hard ball system. 

The base of the induction (i. e. the ergodicity of any two-ball system on a flat 
torus) was proved in [Sin(1970)] and [S-Ch(1987)]. 

Assume now that (M, {5**}, fi) is a given system of (> 3) hard spheres with 
masses mi, m2, . . . , niN and radius r > on the flat unit torus = W jiy (z/ > 2), 
as deflned in §2. Assume further that the ergodicity of every such system is already 
proved to be true for any number of balls AT' with 2 < N' < N. We will carry 
out the induction step by following the strategy for the proof laid down by Sinai 
in [Sin(1979)] and pohshed in the series of papers [K-S-Sz(1989)], [K-S-Sz(1990)-I], 
[K-S-Sz(1991)], and [K-S-Sz(1992)]. 

By using the induction hypothesis, Theorem 5.1 of [Sim(1992)-I], together with 
the slimness of the set A2 of doubly singular phase points, shows that there exists 
a slim subset 5*1 C M of the phase space such that for every x G M \ 5*1 the point 
X has at most one singularity on its entire orbit S^~°°'°°^x, and each branch of 
^(-00, cx))^ is not eventually splitting in any of the time directions. By Corollary 
3.26 and Lemma 4.2 of [Sim(2002)] there exists a locally flnite (hence countable) 
family of codimension-one, smooth, exceptional submanifolds C M such that for 
every point x ^ (Ui ^ '^1 the orbit of x is sufficient (geometrically hyperbolic). 
This means, in particular, that the considered hard ball system (M, {S'*}teR, //) is 
fully hyperbolic. 

By our standing assumption the Chernov-Sinai Ansatz (a global hypothesis of the 
Theorem on Local Ergodicity by Chernov and Sinai, Theorem 5 in [S-Ch(1987)], see 
also Corollary 3.12 in [K-S-Sz(1990)-I] and the main result of [B-Ch-Sz-T(2002)]) is 
true, therefore, by the Theorem on Local Ergodicity, an open neighborhood 3 x 
of any phase point x ^ (|J^ Jj) U Si belongs to a single ergodic component of the 
billiard flow. (Modulo the zero sets, of course.) Therefore, the billiard flow 
has at most countably many, open ergodic components Ci, C2, .... 

Remark. Note that theorem 5.1 of [Sim(1992)-I] (used above) requires the induc- 
tion hypothesis as an assumption. 

Assume that, contrary to the statement of our theorem, the number of ergodic 
components Ci, C2, ... is more than one. The above argument shows that, in 
this case, there exists a codimension-one, smooth (actually analytic) submanifold 
J C M \ dM. separating two different ergodic components Ci and C2, lying on the 
two sides of J. By the Theorem on Local Ergodicity for semi-dispersive billiards, 
no point of J can have a sufficient orbit. (Recall that sufficiency is clearly an open 
property, so the existence of a sufficient point y E J would imply the existence of a 
sufficient point y' & J with a non-singular orbit.) By shrinking J, if necessary, we 
can achieve that the infinitesimal Lyapunov function Q{n) be separated from zero 
on J, where n is a unit normal fleld of J. By replacing J with its time-reversed 
copy 
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-J={{q,v) e M\ {q,-v) G j} , 
if necessary, we can always achieve that Q{n) < —c\ < uniformly across J. 

There could be, however, a little diflBculty in achieving the inequality Q{n) < 
across J. Namely, it may happen that Q{nt) = for every t G M. According to 
(7.2) of [Sim(2003)], the equation Q{nt) = (Vt G M) implies that for the normal 
vector Ut of S^{J) at Xt = S^xq one has nt = (0, Wf) for all t G M and, moreover, 
in the view of (7.5) of [Sim(2003)], — Rw^ is the transformation law at any 
collision xt = {qt, Vt) G 5M. Furthermore, at every collision xt = {qt, Vt) G the 
projected tangent vector ViRw~i~ = Viw'^ lies in the null space of the operator K (see 
also (7.5) in [Sim(2003)]), and this means that wq is a neutral vector for the entire 
trajectory S^y, i. e. Wq G A/q {^S^y). On the other hand, this is impossible for the 
following reason: any tangent vector (5g, 5v) from the space A/q {^S^y^ x Ao {^S^y^ is 
automatically tangent to the exceptional manifold J (as a direct inspection shows), 
thus for any normal vector n = {z, w) G of a separating manifold J one has 

{z, w) e Mo (S^y)^ X ^0 {S^f . 
The membership in this formula is, however, impossible with a nonzero vector 

weAfo {S^y)- 

To make sure that the submanifold J is neatly shaped (i. e. it fulfills (2) of 3.1) 
is an obvious task. Condition (3) of 3.1 clearly holds true. We can achieve (4) as 
follows: Select a base point xq & J with a non-singular and not eventually splitting 
forward orbit S^'^'^^xq. This can be done according to the transversality result 3.3 
above (see also 7.12 in [Sim(2003)]), and by using the fact that the points with an 
eventually splitting forward orbit form a slim set in M (Theorem 5.1 of [Sim(1992)- 
I]), henceforth a set of first category in J. After this, choose a large enough time 
T > so that S'^xq ^ 5M, and the symbolic collision sequence Sq = S (S'^'^lxo) 
is combinatorially rich in the sense of Definition 3.28 of [Sim(2002)]. By further 
shrinking J, if necessary, we can assume that 5^( J) n aM = and is smooth 
on J. Choose a thin, tubular neighborhood Uq of J in M in such a way that be 
still smooth across Uq, and define the set 

(4.1) NS (Uo, So) = jx G Uo\ dimAAo (^S^^'^^x^ > l} 

of not So-sufficient phase points in Uq. Clearly, NS (Uq, So) is a closed, algebraic 
set containing J. We can assume that the selected (generic) base point xq & J 
belongs to the smooth part of the closed algebraic set NS ^Uq, So). This guar- 
antees that actually J = NS (Uq, So), as long as the manifold J and its tubular 
neighborhood Uq are selected small enough, thus achieving property (4) of 3.1. 
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Proof of why property (6) of Definition 3.1 can be assumed. 

We recall that J is a codimension-one, smooth manifold of non-sufficient phase 
points separating two open ergodic components, as described in (0)-(3) at the end 
of §3 of [Sim(2003)]. 

Let P be the subset of J containing all points with non-singular forward orbit 
and recurring to J infinitely many times. 

Lemma 4.2. iJii{P) = 0. 

Proof. Assume that A*i(-P) > 0. Take a suitable Poincare section to make the time 
discrete, and consider the on-to-one first return map T : P — > P of P. According 
to the measure expansion theorem for hypersurfaccs J (with negative infinitesimal 
Lyapunov function Q{n) for their normal field n), proved in [Ch-Sim(2006)], the 
measure jii (T[P)) is strictly larger than //i(P), though T[P) C P. The obtained 
contradiction proves the lemma. □ 

Next, we claim that the above lemma is enough for our purposes to prove (6) of 
3.1. Indeed, the set W d J consisting of all points x e J never again returning to 
J (after leaving it first, of course) has positive //i-measure by Lemma 4.2. Select 
a Lebesgue density base point xq ^ W for W with a non-singular forward orbit, 
and shrink J at the very beginning to such a small size around xq that the relative 
measure of in J be bigger than 1 — 10~^. 

Finally, Main Lemma 3.5 asserts that A ^ 0, contradicting to our earlier state- 
ment that no point of J is sufficient. The obtained contradiction completes the 
inductive step of the proof of the Theorem. □ 

Appendix L The Constants of §3 

In order to make the reading of §3 easier, here we briefly describe the hierarchy 
of the constants used there. 

1. The geometric constant — ci < provides an upper estimate for the infini- 
tesimal Lyapunov function Q{n) of J in (5) of Definition 3.1. It cannot be freely 
chosen in the proof of Main Lemma 3.5. 

2. The constant C2 > is present in the upper measure estimate of Lemma 2 
of [S-Ch(1987)], or Lemma 4.10 in [K-S-Sz(1990)-I]. It cannot be changed in the 
course of the proof of Main Lemma 3.5. 

3. The contraction coefficient < Lq <^ 1 plays a role all over §3. It must 

be chosen suitably small by selecting the time threshold C3 ^ 1 large enough (see 
Corollary 3.20), after having fixed C/q, ^o, and J. The phrase "suitably small" for 
Lq^ means that the inequality 
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^ 0.25//i(J*) 
° 1.2c2C4diam(J) 

should be true, see the end of §3. 

4. The geometric constant C4 > of (3.21) bridges the gap between two distances: 
the distance dist(y, J) between a point y e Dq and J, and the arc length Iq {^p^ ^^^^^ • 
It cannot be freely chosen during the proof of Main Lemma 3.5. 

Appendix II 
Expansion and Contraction Estimates 



For any phase point a; e M \ with a non-singular forward orbit S^^'°°^x 
(and with at least one collision, hence infinitely many collisions on it) we define the 
formal stable subspace E^{x) C T^M of x as 

E\x) = {{5q, 5v) e %M.\ 6v = -B{x)[5q]], 

where the symmetric, positive semi-definite operator B{x) (acting on the tangent 
space of Q at the footpoint where x = {q, v)) is defined by the continued fraction 
expansion introduced by Sinai in [Sin(1979)], see also [Ch(1982)] or (2.4) in [K- 
S-Sz(1990)-I]. It is a well known fact that E^{x) is the tangent space of the local 
stable manifold 7*(a;), if the latter object exists. 

For any phase point a; G M \ 9M with a non-singular backward orbit S'*-~°°''^'*x 
(and with at least one collision on it) the unstable space E'^{x) of x is defined as 
—E^{—x), where —x — (g, —v) for x = {q, v). 

All tangent vectors {5q, 5v) considered in Appendix II are not just arbitrary 
tangent vectors: we restrict the exposition to the "orthogonal section", i. e. to the 
one-codimensional subspaces of the tangent space consisting of all tangent vectors 
{5q, 5v) e XcM. for which 5q is orthogonal to the velocity component v of the phase 
point X = (q, v). 

Theorem. For any phase point xq G M \ dM. with a non-singular forward or- 
bit S^^'°°^xo and with infinitely many consecutive, connected collision graphs on 
S^^'°°^xo, and for any number L > one can find a time t > and a non-zero 
tangent vector {Sqo,Svo) e E^{xo) with 

\\{5qu5v,)\\ 
WiSqo, Svo)\\ 

where {6qt, Svt) = DS\6qo, 6vo) G E^{xt), Xt = S^xq. 
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Proof. Taking time-reversal, we would like to get useful lower estimates for the ex- 
pansion of a tangent vector (Sqo, Svq) G T^qM with positive infinitesimal Lyapunov 
function Q{Sqo, Svq) = {Sqo.Svo). The expression {Sqo,dvo) is the scalar product 
in M'^ defined via the mass (or kinetic energy) metric, see §2. It is also called the 
infinitesimal Lyapunov function associated with the tangent vector {dqo,6vo), see 
[K-B(1994)], or part A.4 of the Appendix in [Ch(1994)], or §7 of [Sim(2003)]. For a 
detailed exposition of the relationship between the quadratic form Q( . ), the rele- 
vant symplectic geometry of the Hamiltonian system and the dynamics, please also 
see [L-W(1995)]. 

Note. The idea to use indefinite quadratic forms to study Anosov systems belongs 
to Lewowicz, and in the case of non-uniformly hyperbolic systems to Wojtkowski. In 
particular, the quadratic form Q was introduced into the subject of semi-dispersive 
biUiards by Wojtkowski, implicitly in [W(1985)], explicitly in [W(1988)]. The sym- 
plectic formulation of the form Q appeared first in [W(1990)]. 

These ideas have been explored in detail and further developed by N. I. Chernov 
and myself in recent personal communications, so that we obtained at least linear 
(but uniform!) expansion rates for submanifolds with negative infinitesimal Lya- 
punov forms for their normal vector. These results are presented in our recent joint 
paper [Ch-Sim(2006)]. Also, closely related to the above said, the following ideas 
(to estimate the expansion rates of tangent vectors from below) are derived from 
the thoughts being published in [Ch-Sim(2006)]. 

Denote by {dqt, dvt) = {DS^){dqo, Svq) the image of the tangent vector {Sqo, dvo) 
under the linearization DS^ of the map S^. (We assume that the base phase point 
xq — for which {Sqo, Svq) e XcqM — has a non-singular forward orbit.) The time- 
evolution (Sqtj^^Svtj^) I— > (Sqt^^Svt^) (ti < t2) on a collision free segment S^^'^'^'^^xq is 
described by the equations 



( A 1 ) 

Correspondingly, the change Q{5qt^,5vt^) i— > Q{Sqt2, Svt2) in the infinitesimal Lya- 
punov function Q{.) on the collision free orbit segment S^^^'^^^xq is 



(A.2) Q{Sqt„Svt,) = Q{Sqt„Svt,) + {t2-ti)\\Svt,\\'', 

thus Q{.) steadily increases between collisions. 

The passage {Sq^,Sv^) i— (dq^^dv^) through a reflection (i. e. when xt = 
S^xo e dM) is given by Lemma 2 of [Sin(1979)] or formula (2) in §3 of [S-Ch(1987)]: 

5qt = R5q7 , 

(A.3) 

5vt = RSv^ + 2 cos (t)RV*KVSq^ , 
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where the operator R : TQ — > TQ is the orthogonal reflection (with respect to 
the mass metric) across the tangent hyperplane Tq^dQ of the boundary dQ at the 
configuration component qt of Xt = {qt,v^), V : (v^)-^ — > Tq^dQ is the t;^~-parallel 
projection of the or tho complement hyperplane (v^)^ onto Tq^dQ, V* : Tg^dQ —>■ 
(v^)-^ is the adjoint of V (i. e. the -parallel projection of Tq^dQ onto (f^")"*", 
where n{qt) is the inner unit normal vector of dQ at qt G dQ), K : Tq^dQ — > Tq^dQ 
is the second fundamental form of the boundary dQ at qt (with respect to the field 
n(g) of inner unit normal vectors of dQ) and, finally, coscj) — {n{qt), vf) > is 
the cosine of the angle (/> (0 < (/> < 7r/2) subtended by vf and n{qt). Regarding 
formulas (A. 3), please see the last displayed formula in §1 of [S-Ch(1987)] or (i)-(ii) 
in Proposition 2.3 of [K-S-Sz(1990)-I]. The instanteneous change in the infinitesimal 
Lyapunov function Q{Sqt,Svt) caused by the reflection at time t is easily derived 
from (A.3): 



Q{Sq+,dv+) = Q{Sq^,dvr) + 2 cos 4>{Vdq^,KVSq^) 
> Q{Sqt ,Svt ). 

In the last inequality wc used the fact that the operator K is positive semi-definite, 
i. e. the billiard is semi-dispersive. 

We are primarily interested in getting useful lower estimates for the expansion 
rate | 1 1/| |5go| I- The needed result is 

Proposition A. 5. Use all the notations above, and assume that 

{5qo, Svo)/\\5qo\\'^ > cq > 0. 
We claim that | |5Qt 1 1/| |5qo| I > 1 + cot for aU t > 0. 



Proof. Clearly, the function \ \Sqt\ \ of t is continuous for alH > and continuously 

A 



differentiable between collisions. According to (A.l), = Svt, so 



(A.6) -\\Sqt\\^ = 2{Sqt,5vt). 

Observe that not only the positive valued function Q{Sqt,Svt) = {Sqt,Svt) is 
nondecr easing in t by (A. 2) and (A. 4), but the quantity {Sqt, Svt) /\\Sqt\ \ is nonde- 
creasing in t, as well. The reason is that {Sqt, Svt) /\\dqt\\ = \\dvt\\cosat {at being 
the acute angle subtended by 6qt and dvt), and between collisions the quantity \ \Svt\\ 
is unchanged, while the acute angle at decreases, according to the time-evolution 
equations (A.l). Finally, we should keep in mind that at a collision the norm \ \Sqt\\ 
does not change, while {5qt,5vt) cannot decrease, see (A. 4). Thus we obtain the 
inequalities 
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{5qt,5vt)/\\8qt\\ > {5qQ,5vo)/\\5qQ\\ > CoH^^oH, 

so 

j^\\5qt\? = 2\\5qt\\j^\\5qt\\ = 2{5qt,Svt) > 2co\\5qo\\ ■ \\Sqt\\ 

by (A. 6) and the previous inequality. This means that J^H^^tH > co||5go||, so 
W^QtW > ||'^<?o||(l + cot), proving the proposition. □ 

Next we need an effective lower estimate cq for the curvature {6qo, 5fo)/||5go|P 
of the trajectory bundle: 

Lemma A. 7. Assume that the perturbation {dq^ , Svq) G T^qM. (as in Proposition 
A. 5) is being performed at time zero right before a collision, say, uo = (1, 2) taking 
place at that time. Select the tangent vector (Sq^ , Svq) in such a specific way that 
SqQ = {m2Wj —miw, 0, 0, . . . , 0) with a nonzero vector w eM.^ , {w, — — 0. 
This scalar product equation is exactly the condition that guarantees that 5qQ 
be orthogonal to the velocity component v~ — (f ^, f^, . . . , f j^) of xq — {q,v~). 
The next, though crucial requirement is that w should be selected from the two- 
dimensional plane spanned by — and qi — q2 (with — q'2|| = 2r) in W^. 
If Vi — V2 and qi — q2 are parallel, then we do not impose this condition, for in 
that case there is no "astigmatism" . The purpose of this condition is to avoid the 
unwanted phenomenon of "astigmatism" in our billiard system, discovered first by 
Bunimovich and Rehacek in [B-R(1997)] and [B-R(1998)]. Later on the phenome- 
non of astigmatism gathered further prominence in the paper [B-Ch-Sz-T(2002)] as 
the main driving mechanism behind the wild non-differentiability of the singularity 
manifolds (at their boundaries) in semi-dispersive billiard systems with a config- 
uration space dimension bigger than 2. Finally, the last requirement is that the 
velocity component 5vq (right before the collision (1,2)) is chosen in such a way 
that the tangent vector (SqQ , Svq) belongs to the unstable space E'^{xo) of a^o- This 
can be done, indeed, by taking Svq = B'^{xo)[SqQ], where B'^{xo) is the curvature 
operator of the unstable manifold of xq at a^o, right before the collision (1,2) taking 
place at time zero. 
We claim that 



+ ) ^ 11^1-^211 ^ \\V1-V2\\ 

W^QoW^ ~ r cos 00 ~ r 

for the post-collision tangent vector (dq^ , 5vq), where cpQ is the acute angle sub- 
tended by — V2 and the outer normal vector of the sphere € M'^l = 2r} 
at the point y — qi — q2- Note that in (A. 8) there is no need to use -|- or — in | |5go| P 
or \\vi -V2\\, for ||5g^|| = \\5q^\\, \\v^ -v^W = \\vt -v^W- 
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Proof. The proof of the first inequality in (A. 8) is a simple, elementary geometric 
argument in the plane spanned by I'l — and qi—q2, so we omit it. We only note 
that the outgoing relative velocity —V2 is obtained from the pre-collision relative 
velocity — V2 by reflecting the latter one across the tangent hyperplane of the 
sphere |y G \\y\ \ = 2r} at the point y = qi — q2- It is a useful advice, though, 
to prove the first inequality of (A. 8) in the case 5vq — first (this is an elementary 
geometry exercise), then observe that this inequality can only be further improved 
when we replace 5vq = with Svq = B'^{xo)[SqQ]. □ 

The previous lemma shows that, in order to get useful lower estimates for the 
"curvature" {Sq,dv)/\\Sq\\'^ of the trajectory bundle, it is necessary (and sufficient) 
to find collisions a = on the orbit of a given point € M with a "relatively 
big" value of \\vi — Vj\\. Finding such collisions will be based upon the following 
result: 

Proposition A. 9. Consider orbit segments S^^''^^xo of A^-ball systems with masses 
mi, m2, . . . , mjv in (or in R^) and with collision sequences E = (cti, cr2, . . . , cr„) 
corresponding to connected collision graphs. (Now the kinetic energy is not neces- 
sarily normalized, and the total momentum X^i^i ^i'^^i may be different from zero.) 
We claim that there exists a positive-valued function /(a; mi, m2, . . . , mAr) (a > 0, 
/ is independent of the orbit segments S^'^'^^Xq) with the following two properties: 

(1) If ||'yi(t/) — Vj{ti)\\ < a for all collisions ct; = (i,j) (1 < ^ < is the time 
of ai) for some trajectory segment S^'^''^^Xo with a symbolic collision sequence E = 
(cTi, a2, . . . , (Tn) corresponding to a connected collision graph, then the norm | \vi' {t) — 
Vj' (t) 1 1 of any relative velocity at any time t e R is at most /(a; mi, . . . , mjv); 

(2) lima^o /(a; mi, ... , mAr) = for any (mi, . . . , ttin)- 

Proof. We begin with 

Lemma A. 10. Consider an AT-ball system with masses mi, . . . , m^v (an (mi, . . . 
, mAr)-system, for short) in T'^ (or in R^). Assume that the inequalities ||vi(0) — 
^-,(0)11 < a hold true {I < i < j < N) for all relative velocities at time zero. We 
claim that 



(A.ll) \\vi{t) -Vj{t)\\<2a 

for any pair and any time t e R, where M = X^^Li "^i 

m = min {mj | 1 < i < N} . 

Note. The estimate (A.ll) is not optimal, however, it will be sufficient for our 
purposes. 
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Proof. The assumed inequalities directly imply (by using a simple convexity ar- 
gument) that 11^^(0)11 < a (1 < i < N) for the velocities v^(0) measured at time 
zero in the baricentric reference system. Therefore, for the total kinetic energy Eq 
(measured in the baricentric system) we get the upper estimation Eq < ^Ma^^ and 
this inequality remains true at any time t. This means that all the inequalities 
||^'i(^)lP ^ ^^"^ hold true for the baricentric velocities vl{t) at any time t, so 

Hit) -v'j{t)\ \ < aVM (m~^^^ + m~^^^^ <2a 
thus the inequalities 

\\vi{t) -Vj{t)\ \ < 2a 

hold true, as well. □ 

Proof of the proposition by induction on the number N. 

For = 1 we can take /(a; mi) = 0, and for = 2 the function /(a; mi, 7712) — a 
is obviously a good choice for /. Let N > 3, and assume that the orbit segment 
xq of an (mi, . . . , mAr)-system fulfills the conditions of the proposition. Let 
(^k = ihj) be the collision in the symbolic sequence T,n — (di, . . . , an) of S^'^''^^xo 
with the property that the collision graph of = (cti, . . . , a^) is connected, while 
the collision graph of Sfc_i = (cji, . . . , a^-i) is still disconnected. Denote the two 
connected components (as vertex sets) of Sfc_i by Ci and C2, so that i e Ci, 
j e C2, Ci U C2 = {1, 2, . . . , iV}, and Ci n C2 = 0. By the induction hypothesis 
and the condition of the proposition, the norm of any relative velocity Vi'{tk — 
0) — Vj'itk — 0) right before the collision ak (taking place at time tk) is at most 
a + /(a; Ci) + /(a; C2), where Ci stands for the collection of the masses of all 
particles in the component C^, Z = 1, 2. Let g{a; mi, . . . , mjv) be the maximum of 
all possible sums 

a + /(a; Di)+f{a-D2), 

taken for all two-class partitions (Di, D2) of the vertex set {1, 2, . . . , N}. According 
to the previous lemma, the function 

fW 

f{a; mi, ... , uin) := 2\ — g{a; mi, ... , uin) 

V m 

fulfills both requirements (1) and (2) of the proposition. □ 

Corollary A. 12. Consider the original (mi, . . . , mAr)-system with the standard 
normalizations X^^Li '^'i'Vi = 0, ^ X^^Li "^iH'^'ilP = |- We claim that there exists a 
threshold G = G{mi, . . . jTUn) > (depending only on N, mi, . . . ,mjv) with the 
following property: 
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In any orbit segment S^^''^^xo of the (mi, . . . , miv)-system with the standard 
normahzations and with a connected coUision graph, one can always find a coUision 
a = taking place at time t, so that — 'I'jl^)!! > G{mi, . . . jTUn)- 

Proof. Indeed, we choose G = G(mi, . . . , mjv) > so small that f{G; mi, ... , mjv) 
< M Assume, contrary to A. 12, that the norm of any relative velocity Vi — Vj 
of any collision of S^'^''^^xo is less than the above selected value of G. By the 
proposition, we have the inequalities ||wi(0) — fj(0)|| < /(G; mi, . . . , mAr) at time 
zero. The normalization Xlili "^i^i(O) — 0; with a simple convexity argument, 
implies that | |f i(0) 1 1 < /(G; mi, . . . , mAr) for alH, I < i < so the total kinetic 
energy is at most |M [f{G; mi, . . . , ttin)]^ < ^, a contradiction. □ 

Corollary A. 13. For any phase point xq with a non-singular backward trajec- 
tory S^~°°'^^xo and with infinitely many consecutive, connected collision graphs on 
S^~°°'^^xo, and for any number L > one can always find a time — t < and a 
non-zero tangent vector {dqo, dvo) G E'^{x-t) {x-t = S~^xo) with | |(5gt 1 1/| |5go| | > L, 
where idqt,Svt) = DS\SqQ,Svo) G ^"(xq). 

Proof. Indeed, select a number t > so big that l-|-^G(mi, . . . ^uin) > L and —t is 
the time of a collision (on the orbit of xq) with the relative velocity v~{—t) —vJ (— t) , 
for which \ \v~ {—t) —v~{—t)\\ > G(mi, . . . , mAr). Indeed, this can be done, thanks 
to A. 12 and the assumed abundance of connected collision graphs. By Lemma A. 7 
we can choose a non-zero tangent vector {Sq^ , Svq) G E'^{x-t) right before the 
collision at time —t in such a way that the lower estimate 

{Sq-^, Sv-^) , 1 s 

holds true for the "curvature" (^Qq", <^f^o~)/||<5Q'o"|P associated with the post-collision 
tangent vector {Sq^ ^ ^'^q)- According to Proposition A. 5, we have the lower esti- 
mate 



||M>i + lG(mi,...,mAr)>L 

OOn r 



|<^9o| 

for the 5q'-expansion rate between {Sq^ , ^Vq") and {5qt, 5vt) = DS*{5qQ, Svq). □ 

We remind the reader that, according to the main result of [B-F-K(1998)], 
there exists a number eo = eo("^i, ■ ■ ■ ,mN', r; v) > and a large threshold A^o = 
No{mi, . . . ,mN] r; u) E N such that in the (mi, . . . ,mN', r; z/)-billiard flow amongst 
any A'^o consecutive collisions one can always flnd two neighboring ones separated 
in time by at least eo- Thus, for a phase point X-t at least eo/2-away from 
collisions, the norms \\5qo\\ and a/I |5qo| P + I l^'f^ol P ai'e equivalent for all vectors 
(Sqo, Svq) G E'^{x-t)i hence, from the previous corollary we immediately get 
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Corollary A. 14. For any phase point xq G M \ with a non-singular back- 
ward trajectory S^~°°'^^xq and with infinitely many consecutive, connected collision 
graphs on S'*-^°°'^-'a;o, and for any number L > one can always find a time —t<0 
and a non-zero tangent vector {Sqo, Svq) e E^{x-t) {x-t = S~^xo) with 



\\{Sqo, Svq] 

where {Sqt,Svt) = DS\6qo,Svo) G £^"(a;o). 



>L, 



□ 



The time-reversal dual of the previous result is immediately obtained by replacing 
the phase point xq = (go, vq) with —xq = (go, — 'i^o), the backward orbit with the 
forward orbit, and the unstable vectors with the stable ones. This dual result is 
exactly our theorem, formulated at the beginning of Appendix II. 
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